Violation of a Bell-like inequality for spin-energy entanglement in neutron polarimetry 
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Violation of a Bell-like inequality for a spin-energy entangled neutron state has been demonstrated 
in a polarimetric experiment. The proposed inequality, in Clauser-Horne-Shimony-Ifolt (CHSIf) 
formalism, relies on correlations between the spin and energy degree of freedom in a single-neutron 
system. The correlation function Schsh is determined to be 2.333 ±0.002 ^ 2, which violates the 
Bell-like CHSH inequality by more than 166 standard deviations. 
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Einstein, Podolsky and Rosen (EPR) argued, based on 
the assumption of local realism, that quantum mechan- 
ics (QM) is not a complete theory In 1951 Bohm 
reformulated the EPR argument for spin observables of 
two spatially separated entangled particles to illuminate 
the essential features of the EPR paradox 0] ■ Thereafter 
Bell introduced inequalities for certain correlations which 
hold for the predictions of any local hidden variable the- 
ory (LHVT) applied Q, but are violated by QM. From 
this, one can conclude that QM cannot be reproduced by 
LHVTs. Five years later Clauser, Horne, Shimony and 
Holt (CHSH) reformulated Bell's inequalities pertinent 
for the first practical test of quantum non locality 
Polarization measurements with correlated photon pairs 
01, produced by atomic cascade [1, 0] and parametric 
down-conversion of lasers [1, d, [l3| , demonstrated viola- 
tion of the CHSH inequality. Up to date many physical 
systems ^ [H, [H llJI have been examined, including 
neutrons in an interferometric experiment [l5| . 

Most EPR experiments test LHVTs, which are based 
on the assumptions of locality and realism. LHVTs are 
a subset of a more general class of hidden-variable the- 
ories, namely the noncontextual hidden- variable theories 
(NCHVTs) . Noncontextuality implies that the value of a 
dynamical variable is determined and independent of the 
experimental context, i.e. of previous or simultaneous 
measurements of a commuting observable Non- 
contextuality is a more stringent demand than locality 
because it requires mutual independence of the results 
for commu ting observables even if there is no spacelike 
separation [T^. 

In the case of neutrons entanglement is not achieved 
between particles but between different degrees of free- 
dom (DOE). Since the observables of one Hilbert spaces 
(HS), describing a certain DOE, commute with observ- 
ables of a different HS, the single neutron system is suit- 
able for studying NCHVTs with multiple DOF. Using 
neutron interferometry [l9| , single-particle entanglement 
between the spinor and the spatial part of the neutron 
wavefunction [l5| , as well as full tomographic state anal- 
ysis [l^l, have already been accomplished. In addition, 
the contextual nature of quantum theory [2l|, [S^l has 
been demonstrated. In a recent experiment creation of 
a triply entangled single neutron state was achieved by 



developing a coherent-manipulation method for the total 
energy (the sum of kinematic and potential energies) of 
a single neutron system [23] . 

Neutron polarimetry has several advantages compared 
to perfect crystal interferometry. It is insensitive to am- 
bient mechanical and thermal disturbances and therefore 
provides better phase stability. Efficiencies of the manip- 
ulations, including state splitting and recombination, are 
considerably high (typically >98 %) resulting in a better 
contrast compared to interferometry. In addition. Single- 
crystal interferometers accept neutrons only within an 
angular range of a few arcseconds, which leads to a sig- 
nificant decrease in intensity. Neutron polarimetry has 
been used to demonstrate the noncommutation proper- 
ties of the Pauli spin operator [2^, and for geometric 
phase measurements [l^, [13, HI] . Polarimetry is applica- 
ble to other quantum s yst ems aside from neutrons, for in- 
stance Ca [29] or '^He [30] ■ Further developments within 
the scope of atomic polarimetry are anticipated. 

This letter reports the first experimental confirmation 
for the violation of a Bell-like CHSH inequality in the 
field of neutron polarimetry. The inequality is based on 
the statistical outcome of correlation measurements of 
the spin and energy DOF. 

In our experiment the neutron's wavefunction exhibits 
entanglement between the spinor and energy DOF [23j . 
expressed as 
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\Ea + nw)®\]) + \E^-kuj)®\i)]. (1) 



Here ]t) and \l) denote the neutron's up and down spin 
eigenstates, referring to the chosen quantization axis, de- 
scribing a spin superposition prepared by a DC-7r/2 spin- 
turner. \Eq + huj) and \Eq — fuj) are the energy eigen- 
states after interaction with a time-dependent magnetic 
field within a radio frequency (RF)-ffipper, due to ab- 
sorption or emission of photons of energy hu!. Eq is the 
initial total energy of the neutron and uj is the frequency 
of the oscillating magnetic field. 

As in common Bell experiments a joint measurement of 
two observables, i.e. spin and energy, is performed: First 
projection operators onto a spin superposition state. 
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FIG. 1: Experimental apparatus for observation of quantum correlations between the spin and energy degree of freedom, 
expressed by a Bell-like ClfSH inequality. The incident neutron beam is polarized by a supermirror polarizer. A DC-n/2 spin- 
turner creates a spin superposition followed by an RF-flipper, preparing the entanglement of spin and energy. The position 
change of the translation stage (displacement of second RF-flipper together with DC spin-flipper) and the phase difference 
between the two oscillating RF-fields adjust the parameters a and 7 for the Bell measurement. Finally, the spin is projected 
back to the initial direction using a second DC-7r/2 spin-turner for a spin polarization analysis followed by a count rate 
detection. The Bloch-sphere description depicts the evolutions of each quantum state in spin and energy subspaces. It includes 
measurement settings of a and 7, determining the projection operators used for joint measurement of spin and energy. The 
effect of the DC-flipper on the translation stage is to suppress a change in the total Larmor precession angle (see main text). 



specified by an angle parameter a, are defined as 
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P'i\a) = _(|t) ±e— |i))((t |±e-(i |). (2) 

Similarly projection operators onto an energy superposi- 
tion state, using an angle parameter 7, are given by 

^ {\Eo + hLu)±e-''^\E^-nw)) (3) 
{{Eo + rvjj\±e''<{Eo~nw\). 



The angle parameters a and 7 are the azimuthal angles 
on the Bloch spheres corresponding to the spin and en- 
ergy DOF, respectively. They are depicted in Fig.[T] In- 
troducing the observables 



and 



i(S)(a)=Pf' (a)-Pi^'(a) 



S(i^)(^) = pf)(^)_p(^^)(^) 



(4) 



(5) 



one can define an expectation value for a joint measure- 
ment of spin and energy along the directions a and 7 

E{a,"f) = (^'(t)|i(^) (g)B(^)|*(i)) =cos(a-f-7).(6) 



For a Bell-like inequality in CHSH-formalism [J four ex- 
pectation values as defined in Eq. ([6]) , with the associated 
directions ai, a2 and 71, 72 for joint measurements of spin 
and energy, are required which yields 



•ScHSif ("1 , ^2 , 71 > 72 ) 
= |£'(ai,7i) -t- E{a2,ji) 



(7) 



E{ai,j2) + £'(02,72) 



The boundary of Eq.([7]) is given by the value 2 for 
any NCHVT, whereas QM predicts a maximal value 
5^^^=2V2 for ai = 0, 02 = n/2, 71 = 7r/4 and 
72 — 37r/4. 

The experiment was carried out at the tangential beam 
port of the 250 kW TRIGA research reactor of the 
Atomic Institute of the Austrian Universities, Vienna. A 
schematic view of the experimental arrangement is shown 
in Fig.[Tl 

A neutron beam of mean wavelength A =1.99 A, re- 
flected from a pyrolytic graphite monochromator and 
propagating in the -|-y-dircction, is polarized along the z- 
direction using a bent CoTi supermirror array. The first 
DC-coil, functioning as a 7r/2 spin-turn device, rotates 
the spin into the xy-plane. Thus a coherent superposi- 
tion of the two orthogonal spin eigenstatesltz) and ||z) 
in equal portions is created, yielding an incident state 
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denoted as 




The entanglement between the spinor and energy DOF 
is created exploiting the mode of operation of a subse- 
quent RF- flipper j23j with an oscillating field B{t) — 

B'f'^ cos{ujt) ■ y. Fulfilhng the resonance condition (lu — 
2\fi\Bo/h) for the oscillating field and the guide field, 
the z-component of the total magnetic field can be com- 
pletely suppressed within the rotating frame of the os- 
cillating field. The effective field, perpendicular to the 
initial polarization, is adjusted to b[^'' = 7r?i/(2r|/i|) ini- 
tiating a spin-fiip process. Here /x is the magnetic mo- 
ment of the neutron and r is the time the neutron re- 
quires to traverse the RF-field region. Interacting with 
a time-dependent magnetic field, the total energy of the 
neutron is no longer conserved due to absorption and 
emission of photons of energy fko, depending on the spin 
state [m, [12 . The RF-flipper is operating at a frequency 
of ti;/27r = 32 kHz and accordingly the guide field is tuned 
to Bq ^ 1.1 mT. The entangled state vector can be rep- 
resented as a Bell state 

I^Bell) = (9) 

-i= {\Eo + fuj) ® IT.) + 1^0 - M ® li.)) , 
which is illustrated using a Bloch sphere description in 

Figm 



The second RF-flipper and an auxiliary DC-flipper are 
mounted on a (single) translation stage. The function 
of the DC-flipper is explained below. By choosing the 
same frequency for the two RF-flippers the energy differ- 
ence between the two spin components is compensated. 
The oscillating field of the second RF-fiipper is given 
by B{t) — b[^^ cos{Lut + (fi^:) • y. This procedure is de- 
scribed by the action of the projection operator for the 
energy recombination = \Eo){{Eo + huj\ + {{EQ - 

fku\). Applied to Eq.®, and considering the second RF- 
flipper(energy recombination) and the DC-flipper, this 
operator yields the final state 

\^nn) = ^ (e-^'^" IT.) + e^^^e'*" U.)) ® \Eo). (10) 

Here w • T = 7 is the phase acquired in energy subspace, 
where T is the propagation time for the distance L + AL 
between the two RF-flippers and is the tuneable phase 
of the oscillating field of the second RF-fiipper. 

The stationary guide field Bq ■ z induces an additional 
phase due to Larmor precession within the guide field 
region. However this phase contribution remains con- 
stant during the experiment and can therefore be ad- 
justed by finding the zero-position of 7 scans (displace- 



ment of translation stage). The corresponding polariza- 
tion vector is given by 

Pfin = ( cos (7 -I- a) , sin (7 4- a) , 0) , (11) 

with the spin phase a = 201^, originating from the phase 
of the oscillating field of the second RF-fiipper (a detailed 
description of the spin phase acquired due to a spinor 
evolutions from the north to the south pole of the Bloch 
sphere, and back, is given in [26|). The phases a and 7 
are associated with the measurement directions on the 
equatorial plane of the Bloch spheres, required for joint 
measurements of spin and energy. 

Compensation of the Larmor phase is accomplished by 
the auxiliary DC-flipper placed subsequently to the sec- 
ond RF-flipper. No additional phase shift, induced by 
Larmor precession, resulting from the change of AL is 
observed. Phase contributions with the same sign occur 
in the regions L + AL (between first and second RF- 
fiippers) and L' — AL (between DC-flipper and second 
DC-7r/2 spin-turner) compensating each other. There- 
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FIG. 2; Typical interference oscillations, due to a variation 
of 7, for ai = 0, a2 — it/2, {oi = tt, = 3tv/2). One 
period corresponds to a displacement of 31.28 ±0.06 mm of 
the translation stage. The dashed lines mark the 7 values 
of 71 = 7r/4, 72 = 37r/4 (7,^ = 57r/4, 7^ = 77r/4 = -7r/4), 
where a maximum violation of the Bell-like CHSH inequality 
is expected. The joint measurements of expectation values 
exhibit Schsh =2.333 ± 0.002. 
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TABLE I: Results of the spin-energy correlation measurements. 







7. (7.^) 


£(Qi,7j) 


(i,i = l,2) 


ai = 0(7r) 


71 


= 7r/4 (57r/4) 


B(ai,7i) 


= 0.594 ± 0.001 


Q2 = 7r/2 (37r/2) 


71 


= 7r/4 (57r/4) 


-E(«2,7i) 


= 0.575 ± 0.001 


ai = 0(7r) 


72 ^ 


= 37r/4 (77r/4) 


-E(ai, 72) 


= -0.571 ± 0.001 


aa = 7r/2 (37r/2) 


72 ^ 


= 37r/4 (77r/4) 


-E(«2,72) 


= 0.593 ± 0.001 



ScHSH= 2.333 ± 0.002 



fore the total Larmor rotation angle ( oc 0Ji^{L + L')) 
remains constant although the position of the translator 
is altered (see [IJl for details of this procedure). Con- 
sequently only 7, the phase of the energy subspace, is 
affected by a displacement of the translation stage. 

The second DC-7r/2 spin-turner reverses the action 
of the first one by a -n/2 spin-rotation around the -fx 
axis for a forthcoming polarization analysis along the 
+z direction by the second supermirror. This is ex- 
pressed by applying a projection operator for the spin 
P^^'' — |Tz)(Tz I- Finally the stationary intensity oscilla- 
tions are given by 

^(«,7) = ^(l + Ccos(a + 7)), (12) 

where C is the contrast, which is 100%(C=1) under 
ideal circumstances. In our setup C was experimentally 
determined as 83.8 ±0.4%, which is depicted in Fig [21 
The physical reasons for the loss in contrast are explained 
later. Thus the expectation value, defined in Eq.®, can 
be rewritten using the normalized count rates obtained 
with the measurement settings of a and 7 denoted as 

_ N{a, 7) + N{a^,j^) - N{a, 7-^) - N{a^,j) 
H", 7) - ^^^^ ^ Ar(a-L, 7^) + N{a, 7-^) -f N{a^,j) ' 

(13) 

with — a + TT and 7^ — •y + n. Therefore, from the 
contrast C=0.838 a value of 2.37(0.838 • 2%/2 - 2.37) is 
expected for Schsh for ai = 0, 02 = 7''/2, 71 = ir/A and 
72 = 37r/4. 

Typical oscillations, observed when the position of the 
translation stage (second RF-flipper) is varied (7-scans), 
are plotted in Fig [2] for different settings of a. One pe- 
riod corresponds to a displacement of the translator stage 
of 31.28±0.06 mm. The 7-scan for ai ~ was used 
to determine the position of the translation stage corre- 
sponding to the values 71 = 7r/4,72 — 3tt/4 (7^*- = 57r/4, 
7^ = 77r/4 = — 7r/4) which are, together with the spin 
phase settings ai = 0, Q!2 = 7r/2 {aj^ = 7r,Q:^ = 37r/2), 
required for determining the S'-value for a maximal vio- 
lation of the Bell-like CHSH inequality. 

The actual Bell measurement consists of successive 
count rate measurements using appropriate settings of 
the phase of the oscillating field of the second RF-flippcr 
(tuning the spin phase a) and position of the second RF- 
flipper mounted on the translation stage (tuning the en- 
ergy phase 7). The four expectation values E(ai,jj) 
= 1,2), for joint measurement of spin and energy 
DOF, are determined from the associated count rates, 



using Eq. 113]). They are listed in Table H] After three 
complete measurement sets (to reduce statistical errors) a 
final value S'cHS_ff=2.333± 0.002 was determined which 
is notedly above the value of 2, predicted by NCHVTs 
and close to 2.37, a value derived by taking a contrast of 
83.8% of the interferograms into account. 

A slight deviation of the measured Schsh value of 
2.333 ±0.002 from the expected value of 2.37 can be ex- 
plained by inhomogeneities of the guide field Bq , which 
results in fiuctuations of the energy phase 7, as well as 
imperfections of the spin phase manipulation. (Vary- 
ing the phase (p^j of the oscillating field of the second 
rf-flipper e.g. by 7r/4 should theoretically yield a phase 
shift of the intensity modulations of 7r/2, as predicted 
by Ea. (|12p . whereas in practice we measured additional 
shifts around one and two degrees from the desired set- 
tings.) Due to the deviation of the measured Schsh 
value from the predicted value it is useful to introduce 
an error estimation of the calculated value of 2.37, which 
consists of three parts: The error of the contrast measure- 
ment (83.8 ±0.4%), inhomogeneities of the guide field 
Bq, and imperfections of the spin phase manipulation. 
These contributions lead to a final error of the calculated 
value estimated by ~ 0.036. This error (2.37 ±0.036) is 
one magnitude larger compared to the error of the Bell 
measurement (2.333 ± 0.002). The former refiects all sys- 
tematic imperfections of the setup, whereas the latter is 
solely a statistic error derived from the count rates, which 
are very high (up to 32000 cuts per point) due to a long 
measurement time. 

The average contrast of C = 83.8%, of the observed 
intensity oscillations, exceeds the minimum visibility 
Ccrit = 70.7%(\/2/2) necessary to exhibit a violation 
of the CHSH inequality. The maximal 5'-value that can 
be achieved experimentally is proportional to the con- 
trast and in our case given by C ■ S^^^g^ — 2.37 > 2 
since C > Ccrit- The main reason for the rather low 
contrast with respect to the high flip efficiencies of the 
DC and RF-flippers is a broad momentum distribution of 
^ 2 % induced by the mosaic structure of our monochro- 
mator crystal. This corresponds to a broad distribution 
of the propagation-time, and therefore the distribution 
of the spin-rotation angle (in the xy-plane) is widened 
after each full rotation of the polarization vector, which 
results a noticeable decrease in contrast. 

In our setup neutrons are detected by a BF3 detector 
with an inherent efficiency of >99 %, which is larger than 
the well known threshold efficiency r/crit = 2(-\/2 — 1) ~ 
0.83 required to close the detection loophole with maxi- 
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mally entangled states [s^, [s^l ■ 

To summarize, we have demonstrated a violation of a 
Bell-like CHSH inequality in the field of polarimetry with 
massive spin- i particles i.e. neutrons. The measured cor- 
relations between the neutron's spinor and energy DOF 
contribute to a result Schsh= 2.333 ±0.002, which ex- 
ceeds the Bell limit, demarcating NCHVTs, by more than 



166 standard deviations. The technique established here 
will be utilized in forthcoming neutron optical experi- 
ments for preparation of multi-entanglement in neutron 
polarimetry and interferometry. 

This work has been partly supported by the Austrian 
Science Foundation, FWF (P21193-N20). 
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